x = -3:.1:3;
px = normpdf(x, 5, 3/2);

xSample = normrnd(5, 3/2, 1e6, 1);
ySample = xSample.^2;

histogram(ySample, 'Normalization', 'pdf', 'BinLimits', [0 80])